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^ ; Abstract 

(^ ' We study four flavor QCD at nonzero temperature and density by analytic continuation from 

^^ ■ an imaginary chemical potential. The explored region is T = 0.95Tc < T < 3.5Tc, and the 

^ . baryochemical potentials range from to ~ 500 MeV. Observables include the number density, 

the order parameter for chiral symmetry, and the pressure, which is calculated via an integral 

method at fixed temperature and quark mass. The simulations are carried out on a 16^ x 4 lattice, 

and the mass dependence of the results is estimated by exploiting the Maxwell relations. In the 



(N 



^Sj ■ hadronic region we confirm that the results are consistent with a simple resonance hadron gas 

^ . model, and we estimate the critical density by combining the results for the number density with 

,__! i those for the critical line. In the hot phase, above the endpoint of the Roberge- Weiss transition 

^D I Te — 1.17c the results are consistent with a free lattice model with a fixed effective number of flavor 

^— N ' slightly different from four. We confirm that confinement and chiral symmetry are coincident by a 

■^ ■ further analysis of the critical line, and we discuss the interrelation between thermodynamics and 



o 



critical behavior. We comment on the strength and weakness of the method, and propose further 



C^ , developments. 
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I. INTRODUCTION 

QCD at nonzero temperature and density is an important subject both from a theoretical 
and phenomenological perspective y|. Current and future experiments at RHIC and LHC 
will explore the region of the phase diagram close to the zero density axis, not far from 
the cooling path of the primordial universe. The cold, high density phases are relevant for 
astrophysics. Future GSI experiments might bridge these two regimes and explore the region 
of intermediate temperatures and densities. 

While many predictions on the structure of the phase diagram can be obtained by use of 
simple models dictated by symmetry arguments, quantitative studies require a first principle 
calculation. Lattice field theory is the natural approach, and QCD at nonzero temperature 
and density is by now studied by a variety of lattice methods: the multiparameter reweight- 
ing method achieves on optimal overlap between the simulation ensemble at zero baryon 
density and the target ensemble at nonzero density [2, Isl y, |5| ; the direct calculation of the 
derivatives gave the first informations on the physics at nonzero chemical potential p, LD, |8[ ; 
the Taylor expanded reweighting reduces the numerical costs associated with the multipa- 
rameter reweighting |3, ilO] ; the analytic continuation from an imaginary chemical potential 
uses the information from the /^^ < halfplane to reconstruct the physics of real baryon 
density ll|, |l2, llM [ij, ll5| . In addition to these pragmatic works, which take advantage of 



the fluctuations close to Tc to explore the small chemical potential region fi/T < 1., there 
have been a number of new proposals (I£,il7|,|l8|, discussions |.19j a nd checks |20]. For recent 
reviews see 2l|,|2j|, and for introductions into the subject see J23|, |2J. 

In this paper we extend our study of four flavor QCD within the imaginary chemical 
potential approach. Let us remind that four flavor QCD has a flrst order transition at 
/i = 0, as it is the one at T = 0. Hence, it is expected to have a flrst order critical line 
in the plane T, /i: there is no tricritical point or endpoint to be investigated here. On the 
other hand, the flrst order (or sharp crossover) nature of the critical line makes the model 
simpler, and amenable to a detailed study with comparatively modest numerical resources. 

In our flrst paper \l4\ we studied the order parameter in the hadronic phase and calculated 
the critical line up to /x ~ 500MeV. We have indeed conflrmed the expected flrst order 
(or very sharp crossover) nature of the transition, we have observed its interrelation with 
deconflnement, and we have studied its Nf scaling by comparing our results with those of ref. 
[l3| . The results on the chiral condensate were consistent with a partition function described 
by a simple hyperbolic cosine behavior. Consistent numerical flndings were reported in ref. 
J31[, and intepreted within a resonance gas model. 

This work addresses more fully the properties of the hadronic phase, and those of the 
plasma phase. We base our analysis mostly on the results for the number density and 
for the chiral condensate. We calculate the subtracted pressure by an integral method 
at constant temperature and quark mass, and we study the mass dependence by taking 
numerical derivatives of the chiral condensate. 

The properties of the quark gluon plasma phase are studied by means of the number 
density and of the subtracted pressure. We discuss the nonperturbative nature of the hot 
phase in the vicinity of the critical point, and we argue that this result follows naturally by 
an analysis of the phase diagram in the T, /z^ plane. We compare the numerical results with 
analytic predictions, a task for which the imaginary chemical potential approach is ideally 
suited: for instance, rather than analytically continue the numerical results to real chemical 
potential, we can continue the analytic predictions to the negative /x^ halfplane, and contrast 



the results with the numerical ones. 

In the hadronic phase we study in detail the number density: its simple behavior will 
further support the applicability of the hadron resonance gas model; moreover we combine 
the results for the critical line with those for the number density in order to estimate the 
critical density. 

In both phases we exploit the Maxwell relations to study the mass dependence: this turns 
out to be sizable in the hadronic phase, and negligible in the plasma phase. 

Concerning the critical line, as already mentioned, the four flavor model has a first order 
transition. Hence, tricritical points or endpoints are not present here, and we contented 
ourselves with the precision reached in our previous study. The new results on the crit- 
ical behavior consists in a more detailed analysis of the correlation of the chiral and the 
deconfinement transition. 

The paper is organised as follows. In Section II we review the properties of the phase 
diagram of QCD in the chemical potential- temperature plane, including the new results 
on the critical behavior at a selected fi value. In Section III we describe our observables 
and give an overview of the results. Section IV is devoted to the results in the hadronic 
region, while Section V presents results for the Quark Gluon Plasma phase. In either phases 
we will discuss in detail the various ansdtze which emerge naturally once the analyticity 
properties and the nature of the critical lines are taken into account. In Section VI we 
discuss our results for Tc < T < Te and give some general comment about this intermediate 
region of temperatures. We summarize and discuss future perspectives in Section VII. Some 
preliminary results have already appeared in [25[. 



II. QCD IN THE /u^ T PLANE 

Results from simulations with an imaginary chemical potential can be analytically con- 
tinued to a real chemical potential, thus circumventing the sign problem [Ij, |l3|, |lj|. In 
practice, the analytical continuation is carried out along one line in the complex /x plane: 
first along the imaginary axes, and then along the real one. It is then meaningful to map 
this path in the complex /i^ plane: because of the symmetry property Z{fi) = Z{—ji) this 
can be achieved without losing generality. In the complex /x^ plane the partition function is 
real for real values of the external parameter /i^, complex otherwise: the situation resembles 
that of ordinary statistical models in an external field. Hence, the analyticity of the physical 
observables [ll| as well as that of the critical line [l3| follows naturally. 

The reality region for the partition function represents states which are physically acces- 
sible. The reality region for the determinant represents the region which is amenable to an 
importance sampling calculation: 3fJe/i^ < 0. 

The phase diagram in the temperature, (real) yU^ plane is sketched in Fig.l, where we omit 
the superconducting and the color flavor locked phase, which (unfortunately) play no role 
in our discussion. The region accessible to numerical simulations is the one with /i^ < 0: at 
a variance with other approaches to finite density QCD, which so far only used information 
at /i = [2, llOlllSl, the imaginary chemical potential method exploits the entire halfspace. 



Note that after the rotation to /i^ there is no analytic continuation in complex plane to 
be done, but rather we can talk about a simpler analytic extrapolation along the real /i^ 
axis. 

We note here that there are physical questions which can be addressed quite simply: in 
Fig. 121 we demonstrate the correlation between the average values of the chiral condensate 




FIG. 1: Sketch of the phase diagram in the /i^,T plane: the sohd line is the chiral transition, the 
dashed line is the Roberge Weiss transition. Simulations can be carried out at ^^ < and results 
continued to the physical domain fi'^ > 0. The derivative and reweighting methods have been 
used so far to extract informations from simulations performed at fj, = 0. The imaginary chemical 
potential approaches uses results on the left hand half plane. Different methods could be combined 
to improve the overall performance 



and of the Polyakov loop. In Fig. El the same correlation is illustrated directly on the 
Monte Carlo time histories at the phase transition, showing the striking result that not only 
the two phase transitions are coincident, but that the chiral condensate ant the Polyakov 
loop are completely correlated even at the level of small scale fluctuations. The data are 
obtained at fixed lit = 0.15, and a similar behavior can be observed for other values of /x/, 
including zero 26||. This correlation should be continued at real baryon density: to this 
end, we note that if /3c(i/i/) = Pd{ifJ'i) over a finite imaginary chemical potential interval, 
then the function Aj3{ifj,i) = Pd'ifJ'i) — l^diil^i) is simply continued to be zero over the 
entire analyticity domain, thus demonstrating the correlation between the chiral and the 
deconfinement transitions (/5c = /3d) also for real values of ji. We refer to [27[ for an effective 
Lagrangian discussion of this issue, and to ^] for analogous results in the two color model. 
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FIG. 2: Correlation between (ipip) and Polyakov loop at /U/ = 0.15, demonstrating the correlation 
of chiral and deconfining transition at nonzero baryon density. 

III. OBSERVABLES AND OVERVIEW OF THE RESULTS 

The numerical simulations were performed, using the HMC algorithm, with four flavor 
of staggered fermions, on the same lattice 16^ x 4 and with the same mass am = 0.05 as 
in our previous study. New numerical results have been obtained in the plasma phase, for 
j3 = 5.310, 5.650, 5.689, corresponding to T ~ (1.5Tc, 2.5Tc, 3.5Tc) (we used the two loop 
P function to convert to physical units), and statistics for (3 = 5.030(T ~ 0.985Tc) and 
P = 5.100(T ~ 1.095Tc), close to the Roberge Weiss transition point, have been improved 
as well. 

Our analysis is mostly based on measurements of the chiral condensate, number density 
and Polyakov loop. Note that the number density is purely imaginary for imaginary chem- 
ical potential: in the following n{T, fii,mq) will denote the imaginary part of the result. 
From (tptp) and n we will calculate and exploit derivatives and integrals with respect to the 
chemical potential, yielding the mass dependence of the number density, the quark number 
susceptibility , and the pressure. 

In Figure |3] we present an overview of our results for the number density, where we can 
read off the main features outlined in the discussion of the phase diagram of Section II above: 
for T < Tc the results are smooth; in the intermediate region Tc < T < Te^, there is a clear 



^ We will use Tg to denote the cndpoint of the Roberge Weiss transition; as there is no QCD endpoint in 
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FIG. 3: Monte Carlo histories of the chiral condensate (upper history) and of twice the Polyakov 
loop modulus (lower history) for /i/ = 0.15 at the phase transition. The apparent correlation 
between the two quantities is the clearest demonstration of the coincidence of deconfinement and 
chiral phase transition also for /i/ ^ 0. 

discontinuity, in correspondence with the chiral/deconfining transition; finally, above Te , 
the results for the number density increase sharply, eventually approaching the free field 
results. 

Slightly anticipating the numerical analysis which will be presented below, we show in 
Figure |31 the summary of our results for the quark number susceptibility Xg'- 



Xq{T,fi = 0,mq) 



dn{T,fii,m g)^ 



Im=o 



(1) 



Obviously, in the case of a polynomial behavior n(T, /x/, rriq) = a(T, mq)fii + b{T, mq)fi^j, 
XqiT,fi = 0,mq) = a(T,mq). For other fits Xq{T,fi = 0,mq) is obtained from equation (jT)), 
taking for n{T, fij,mq) the corresponding fitted function. The agreement between different 
analysis can be judged from Figure El 

One word of comment might be in order: aside from the free massless case, the analytic 
form for n{T, fij, rriq) is, of course, not known. Given that the numerical results have errors. 



this model no confusion should arise 



different ansiitze might well produce equally satisfactory results (measured, for instance, 
by the x^/d.o.f). When this is the case, we have checked that the analytically continued 
results from the two different forms are consistent with each other, within errors. The results 
presented in Figure offer one first example of this check. 

To calculate the pressure, or rather, its variation with /i ^ at constant temperature and 
quark mass 

AP/T^ = {P{T, /i, m,) - P{T, /i = 0, m,))/T^ (2) 

we exploited the relationship 

n{T, /i, rrig) = — (3) 

which gives 

{P{T, /i, m,) - P{T, /i = 0, m,))/T' = Nf J n{fi)dfi (4) 

This can be achieved by numerical integration of the results for the number density 
n(T, fii,mg)), thus obtaining AP(T, /i/,mg)/T^ (which can be continued to real chemical 
potential). A similar approach to the calculation of pressure can be pursued at zero baryon 
number |32|. 

In Figure ini we show the results of AP(T, /^/,mg)/T^ obtained in this way. The results 
for the free case have been obtained by calculating tt, on a 12^ x 4 lattice, for rriq = .05, 
fitting the resulting data to a third order polynomial, then continued to imaginary chemical 
potential, which amounts to a flip of the third order term. Actually, for T > I.IT^ a linear 
term suffices to describe the data, the third order term can be set to zero, and the data 
in Figure IHl coincide with their analytic continuation to a real chemical potential in that 
temperature range. 

In Figure [7| we plot the same data in a different form, for the sake of an easier comparison 
with results from other groups jsl |lO|, |33|, |3^ . We also note that an alternative procedure to 



obtain AP(T, /x, mq)/T'^ is of course by analytical integration of the results for n(yu) obtained 
by analytic continuation: the two procedures give consistent results. 

In J25|, we computed the fermionic contribution to the pressure using the basic thermo- 
dynamic identity 

p(r.^.,n,)^ ™"^'^;J-^-"''> (5) 

combined with the naive tree level result for the Karsch coefficients |29|]. By comparing the 
results for AP/T^ of J25| to our present results from the integral method, we note that the 
nonperturbative Karsch factors should be fairly large: indeed, the tree level results exceed 
the numerical results from the integral method; even ignoring the gluonic contribution, we 
see that the correcting factor should be ~ .5. 

One further observable we shall consider is the derivative of the chiral condensate with 
respect to the chemical potential. This can be computed by numerical differentiation and will 
be used to estimate the mass dependence of the number density according to the Maxwell 
relation bc 



d{'ip'ip){T,ij,,mq) _ dn{T,^i,mq) 



dn dm, 



(6) 
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^ From now on /i will denote the real chemical potential, i.e. the complex chemical potential is ^icompie 
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FIG. 4: Overview of the results for the baryon density as a function of fij: the behavior is smooth 
in the hadronic phase, shows the expected discontinuity associated with the chiral/deconfining 
transition in the intermediate region Tc < T < Te, and increases rapidly in the quark gluon 
plasma phase. 



IV. THE HADRONIC PHASE AND THE HADRON RESONANCE GAS MODEL 

In this region observables are a continuous and periodic function of fii/T , analytic 
continuation in the /i^ > half plane is always possible, but interesting only when 

Xg{T,fl = 0)>0. 

The analytic continuation of any observable O is valid within the analyticity domain, i.e. 
till II < ficiT), where fidT) has to be measured independently. The value of the analytic 
continuation of O at fic, 0{fic), defines the discontinuity at the critical point, or, equivalently, 
its critical value. This allows the identification of the order of the phase transition: first, 
when 0{fic) 7^ 0, second, when 0{fic) = 0. 

Taylor expansion and Fourier decomposition are among the natural parameterizations 
for the observables 1J| . In particular, the analysis of the phase diagram in the temperature- 
imaginary chemical potential plane suggests to use Fourier analysis in this region, as ob- 
servables are periodic and continuous there. Note that for the simple one dimensional QCD 
model, which can be analytically solved and is related to the partition function of QCD in 
the infinite coupling limit, all of the Fourier coefficients but the first ones will be zero. 

For observables which are even (Og) or odd (Oq) under /i -^ —^ the Fourier series read: 



Oe{fii,Nt) = 5]aScos(niV,iV,/i,) 

n 



(7) 
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FIG. 5: Particle number susceptibility at ;U = as a function of temperature from different fits (see 
text for details). The susceptibility is small, yet different from zero, in the hadronic phase close to 
Tc, and reaches a nearly constant value already for T ~ 1.5Tc. 



n 

which is easily continued to real chemical potential: 

Oe{fii,Nt) = J2(^Flcosh{nNtNcfii) 

n 

Oo{iJ,i,Nt) = ^aSsinh(niV,iV,/i,). 



(8) 

(9) 
(10) 



In our past Fourier analysis of the chiral condensate - obviously an even observable - we 
limited ourselves to n = 0,1,2 and we assessed the validity of the fits via both the value of 
the x^/d.o.f . and the stability of apl and apl given by one and two cosine fits: we found that 
one cosine fit is actually enough to describe the data up to T ~ Tc in the four flavor model 



14| : adding a term cos{2NtNciJ.) = cos(24yu) in the expansion did not modify the value of 



the first coefficients and does not particularly improve the x^/d.o.f.. 
We present here the analogous analysis for the number density. 
A first round of fits was performed by setting all of the Fourier coefficients but the first 

(1) + 

one ap^ = ap to zero 



n(T, fii,mq) = apsin{?)jii/T) 



obtaining 



n{T,ij,j,mg) 
n{T,fii,mq) 



0.0039{2)sin{3fii/T) 
0.0062(2)sm(3/X7/T) 



(11) 

(12) 
(13) 
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FIG. 6: Overview of the results for AP(T, jij, mq)/T'^ as a function of iJ,j from the integral method, 
and different temperatures. 

for (3 = 5.010 (upper, xVd.o.f. = 0.39) and f3 = 5.030 ( xVd.o.f. = 1.06), or equivalently 
T ~ 0.95Tc and T ~ 0.985Tc. The results, and the relative errorbands, are shown in Figure 

El 

To further assess the validity of this simple parametrization, we have also performed 
polynomial fits of the form 

f{x) = Ua^p'^x - 288ag^x^ (14) 



The constants are chosen so that ap 



(3 = 5.010 we find a 
with a 



(1) 



ap = CL 

(3) 



(3) 

p when a simple sine fit is adequate. At 



(1) 



ap 



,(3) 



ip , 



-0.0044(7) and a^ 



-0.009(4) (xVd-o.f. = 0.28) 



within the large errors. At /? = 5.030 we find a 



(1) 



consistent 
-0.0060(2) 



and ap' = —0.0048(3) (x^/d.o.f. = 0.69), again indicating that the the contribution from a 
second Fourier coefficient is small, although perhaps not negligible. 

The obvious analytic continuation, representing the results for the number density for 
these parameters from one parameter Fourier fit, read: 



n{p = 5.010(T ^ 0.95Te), /i, m^) 
n{p = 5.030(T ~ 0.985Tc), /i, rUg) 



0.0039(2)sm/i(3/i/T) 
0.0062(2)sm/i(3/i/T) 



(15) 
(16) 



A second order term in the Fourier series improves slightly the quality of the fit but is 
poorly determined. It is however important to determine the uncertainty on the analytic 
continuation and on the estimate of the critical density induced by this contribution: indeed, 
it might well be that a term which is subleading in the imaginary chemical potential domain 
becomes leading when the results are analytically continued to the real plane. 
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FIG. 7: Overview of the results for AP(T, ni,mq)/T^ as a function of T/Tc, and different values 
of m/T. 



The situation is better seen in the plot Figure El where we plot the results for the analytic 
continuation to real /i with the errorbands from the fits, for one (dotted line) and two 
(dashed) Fourier coefficients fit, for /3 = 5.030. On the same plot we mark with a vertical 
line at /i = /ic the limit of validity of the analytic continuation, which coincides with the 
limit of the hadronic phase. 

Following the discussion above, at /x = /ic we read off the plot n{^c) = 0.0087(30) (lattice 
units). If we consider the second Fourier's coefficients the induced uncertainty grows big, and 
we would obtain n{fic) = 0.011(1). All in all, the critical density we estimate at T = 0.985Tc 
is n{fi)/T^ ~ 0.6. 

To study (at least semiquantitatively) the mass dependence of the results, we consider 
the Maxwell relation 

d{ijjijj)/dfi = dn{^)/dm . (17) 

The results for the chiral condensate can thus be used to estimate the mass dependence. 

In an attempt to introduce as less prejudice as possible, we have first numerically derived 
the results for the chiral condensate. The results, although very noisy, are in agreement with 
the derivative of the fitting function for the chiral condensate itself. So we use the latter for 
the subsequent discussion. 

Consider the parametrization for the chiral condensate 



which combined with 



(V'-?/')(T, /i,mg) = accosh{3fiNT) + be 



n{T, fi,mq) = anSinh{3fiNT) 



(18) 
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FIG. 8: One Fourier coefficient fit to the particle number in the hadronic phase. 



gives 



n(/x, niq + Amg) — n{jj,^ m^ 



SNtAm 



ac 



(20) 



Combining the present results form the number density with the old ones [IJ] for the 
chiral condensate, we obtain 



n{^, Trig + Anig) — n{fi, rrig) 

n{fi, nig) 
n{iJi, Trig + Amg) - n(yU, m^) 

n{iJ,mg) 



2.5 X 3NtAm„ 



4 X 3N.Am„ 



(21) 
(22) 



for T = 0.95Tc {f3 = 5.010) and T = 0.985Te (/3 = 5.030) respectively. 

To illustrate this dependence, we plot the numerical results in Figure ^| Note that in 
the plots fi range from 0.0 to [idT): hence the critical densities and their mass dependences 
can be read off at the intercept with the righthand border of the plot. 

We meant here to get an estimate of the mass dependence of the number density: to this 
end we used for both temperatures only the results from one Fourier coefficient fit, and we 
omit the errors (which, for T = 0.985Tc, can be read off Figure |H1) 

To summaries our findings, our previous results for the chiral condensate jlj] and the 
present ones for the number density are consistent with AP oc {cosh{fiB/T) — 1) in the 
broken phase. There is room for a small deviation especially at higher T values. 

In J31[ it was shown that the data obtained from an expanded reweighting behave in 
the same way, and it was pointed out that this result is the one expected from an hadron 
resonance gas model. 
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FIG. 10: The mass dependence of the number density at T = .985Tc (left) and T = 0.95Tc (right) 

V. THE QGP PHASE AND THE EQUATION OF STATE 

At high temperature, in the weak coupling regime, perturbation theory might serve as a 
guidance, suggesting that the first few terms of the Taylor expansion might be adequate in 
a wider range of chemical potentials. So, at a variance with the expansion in the hadronic 
phase, where the natural parametrization is given by a Fourier analysis, in this phase the 
natural parametrization for the grand partition function is a polynomial. 

The leading order result for the pressurep(T, yu) in the massless limit is easily computed, 
given that at zero coupling the massless theory reduces to a non-interacting gas of quarks 
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TABLE I: Coefficients of a polynomial fit for the number density in the Quark Gluon Plasma 
phase. 

I < 4_ xVd.o.f. 

5.10 -0.4646(68) 2.02 (60) 0.89 

5.310 -0.4994(40) 1.83(64) 0.92 

5.650 -0.5129 (43) 2.36(82) 1.65 

5.869 -0.5087 (16) 0.89(28) 0.30 

5.10 -0.4442(42) 1.66 

5.310 -0.4897(29) 1.87 

5.650 -0.5026 (31) 2.88 

5.869 -0.5039 (7) 0.58 



and gluons, yielding for the pressure 



^2 



^ Ti O , ■,M'"f\ , ""■!.. 2t2 , "/ 



PP-, ^) = ^r' (^8 + 7iV.:j j + '^^^T^ + j^_^^ . (23) 

Obviously, when analytically continued to the negative /i^ side, this gives 



,2 



P(T, ,,) = |r' (8 + 7nJ!±) - '^,]T' + ^y, . (24) 

Because of the Roberge Weiss periodicity this polynomial behavior should be cut at the 
Roberge Weiss transition /i/ = ttT/3: this is consistent with the Roberge Weiss critical line 
being strongly first order at high temperature. We discuss first the results of the fits of the 
number density to polynomial form; then we contrast these results with a free field behavior. 
The considerations above suggests a natural ansatz for the behavior of the number density 
in this phase as a simple polynomial with only odd powers. We performed then fits to 

niT,fxi) = a{T)fxi-b{T)fi'j (25) 

whose obvious analytic continuation is 

n{T, /i) = a{T)i^ + ^(T)/!^ . (26) 

Note again that a(T) = Xq(T,fj, = 0). 

The results of the fits are given in Table |lj upper rows. To assess the relevance of the 
third order term we have performed fits with b{T) = 0, whose results are summarized in 
the bottom rows of Table HI As usual, the quality of the fit worsten slightly, while the first 
coefficient a{T) remains compatible with that estimated by a two parameter fit. 

On the other hand, at /i/ = 0.1 (for instance) the contribution of the third order term to 
the number density of the free lattice gas is below two percent : a fair set of measurements 
with this precision around fij = 0.1 would then be needed to disentangle the third order 
term from the error, and it comes as no surprise that within the current precision is not 
possible to safely estimate it. 

In Figure ^2 we show the results for the particle number at T = 1.5Tc, T = 2.5Tc, 
T = 3.5Tc as a function of the imaginary chemical potential, together with the free lattice 
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result (because of the known discrepancies between the lattice and continuum behavior in 
the free case at Nt = 4, we used lattice free results for this comparison, as was already done 
for the quark number susceptibility in Figure 4 above). 

Some deviations are apparent, whose origin we would like to understand. It would be 
however arduous, given the strong lattice artifacts, to tr y to make contact with a rigorous 
perturbative analysis carried out in the continuum jSy, |37|, |38|. Rather then attempting 
that, we parametrize the deviation from a free field behavior as 3311351 

AP(T,/i) = /(T,/.)P/_(T,/x) (27) 

where Ph^^ (T, p,) is the lattice free result for the pressure. For instance, in the discussion of 
Ref. 

/(T,/i) = 2(1 -2a,/7r) (28) 

and the crucial point was that a^ is fi dependent. 

We can search for such a non trivial prefactor f(T,fi) by taking the ratio between the 
numerical data and the lattice free field result njj.^^{fij) at imaginary chemical potential: 

fi(r.„) = 4^ (29) 

A non-trivial (i.e. not a constant) R{T,fij) would indicate a non-trivial f{T,fi). 

In Fig. [T21we plot R{T, fij) versus fii/T: we see that R{T, fii) is constant within errors, so 
that our data do not permit to distinguish a non trivial factor within the error bars: rather, 
the results for T > 1.5Tc seem consistent with a free lattice gas, with an fixed effective 
number of flavors iV/^(T)/4 = i?(T): Nf^ = 0.92 x 4 for T = 3.5Te, and Nf^ = 0.89 x 4 
for T = 1.5Tc. These results confirm those obtained at /i = from the quark number 
susceptibility in Figure 4 and extend them over a finite range of chemical potentials 

One last remark concerns the mass dependence of the results, which, as in the broken 
phase, can be computed from the derivative of the chiral condensate. In the chiral limit 
this gives ^ = , since the chiral condensate is identically zero. We have verified that 
^ remains very small compared to n itself: in a nutshell, in the quark gluon plasma phase 
< tpip > is very small (zero in the chiral limit), while the number density grows larger, and 
this implies that the mass sensitivity is greatly reduced with respect to that in the broken 
phase. 



VI. THE INTERMEDIATE REGIME T^ < T < Te 

The discussions presented above bring very naturally to the consideration of a dynamical 
region which is comprised between the deconfinement transition, and the endpoint of the 
Roberge Weiss transition. 

In this dynamical region the analytic continuation is valid till fi = oo but the interval 
accessible to the simulations at imaginary fi is small, as simulations in this area hits the 
chiral critical line for /i^ < 0. 

In Figure El we repeat the same analysis for the number density done in the previous 
Section, but for T/Tc = 1.095. We see that the results are noisier than those as higher tem- 
perature, and it is difficult to draw firm conclusions. Anyway, they might still accommodate 
some deviation from a simple free field with a reduced effective number of flavor. 
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FIG. 11: Results for the particle number as a function of imaginary chemical potential. The dotted 
line is the free field for imaginary chemical potential. 



Let us make some general consideration about the thermodynamic behavior in this region 
by considering the critical line at imaginary chemical potential, Let us consider first the 
case of a second order transition: the analytic continuation of the polynomial predicted by 
perturbation theory for positive /i^ would hardly reproduce the correct critical behavior at 
the second order phase transition for /i^ < 0. In fact, for a second order chiral transition 
at negative /i^, AP(T, /i^) oc (jU^ — /U^^)'^, where x is a generic exponent. As the window 
between the critical line and the /i = axis is anyway small, such behavior - possibly with 
subcritical corrections - should persist in the proximity of the real axis. For generic values of 
the exponent a second order chiral transition seems incompatible with a free field behavior. 
The same discussion can be repeated for a first order transition of finite strength, by trading 
the critical point fic with the spinodal point fi* . So deviations from free field are to be 
expected in this intermediate regime. 



VII. SUMMARY AND OUTLOOK 

We have gained a good understanding of the strength and weakness of the method. 
First, the imaginary chemical potential approach is not limited to small volumes (aside 
from the usual limitations of any lattice calculations). Next, physical observables can be 
directly computed by usual methods, and their analytic continuation, or, extrapolation, 
can be pushed up to the critical line, thus providing estimates of the critical values and 
discontinuities. In addition to that, the method provides a natural test bed for analytic 
models or calculations, which can be analytically continued to imaginary chemical potential. 
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the ratio of the lattice results to the lattice free field. 
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and directly contrasted with the numerical results. On the weak side, corrections which are 
subleading for imaginary chemical potential, or for /i^ < 0, might become leading in the real 
domain, /i^ > 0.: it is important to try and cross check different analytic parameterizations, 
and we have given a few examples of this procedure in our analysis. 

We have obtained results on the four flavor model for .985Tc < T < 3.5Tc, and fiB < 
500MeV. 

Concerning the critical line, we have studied in detail the chiral and "deconflning" tran- 
sition at a selected value of ni and conflrmed that they remain correlated, showing also the 
complete correlation between the Monte Carlo time histories of the Polyakov loop and the 
chiral condensate around the phase transition. As explained in Ref. IJ] and in Section II 
above, this, together with the observation of their correlation at zero chemical potential, 
implies the equality of the critical temperature 



T^il^) = Tt{^^) 



(30) 



also for real chemical potential. 

In the hadronic phase the corrections to n(T,fiB) oc sinh{^B/T) are very small. This 
confirms and completes the finding of Ref. |lj| where we did show that the chiral condensate 
behaves as {ipip)(T, fis) oc cosh{^B/T) + c. In conclusion our results in the hadronic phase 
are consistent with an hadron resonance gas model, possibly with small corrections close 
to Tc. Again in the hadronic phase we have calculated the baryon density in the hadronic 
phase, and estimated its critical value n{fic,T = .9S5Tc,mq = .05)/Tj^ ~ 0.6; the mass 
dependence has been inferred from the Maxwell relation giving An = —A.03Amq/T. 

In the high temperature regime, for T > 1.5Tc the results are compatible with lattice 
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FIG. 13: Deviation from the free field behavior in the RW regime 

Stefan-Boltzmann with an effective fixed number of active flavors ~ 0.92 x 4 for T=3.5 T^. 
and ~ 0.89 x 4 for T = 2.5Tc. We found that the mass dependence is very small in this 
region. 

We discussed the interplay between thermodynamics and chiral transition in the region 
comprised between the critical point Tc and the endpoint of the Roberge Weiss transition Tg. 
We noted the possibility of non-trivial deviations from a free lattice field, possibly connected 
with the chiral transition at /i^ < 0. 

As for future applications, the method seems ideally suited for more detailed comparisons 
with analytic models, and , more important, nothing prevents its extension to larger lattices. 

We think that the performance could be further improved by considering hybrid methods 
which combines the imaginary chemical potential approach with other methods, for instance 
by making use of reweighting J3|, |l8| or direct calculations of derivatives p at nonzero fx to 
improve the accuracy of the results at negative /x^. 

Finally, the study of discontinuities, as sketched in Sect. IV above, might offer an alter- 
native approach to the study of endpoints and tricritical points. 
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